Comment on "A Kac-potential treatment of nonintegrable interactions" by 

Vollmayr-Lee and Luijten 
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We com ment the recent m anuscript by Vollmayr-Lee 
and Luijten |cond-mat /0009031 1 focusing on classical systems 
with nonintegrable interactions. The authors claim that they 
have proved that Boltzmann-Gibbs statistics suffices for de- 
scribing the system in thermal equilibrium. We show that this 
statement is a misleading oversimplification since it only ap- 
plies for the limjv^cx) limt_ 00 ordering, but certainly not for 
the limt^oo limjv^oo one. The latter can even be the unique 
physically meaningful situation for thermodynamically large 
systems. 

In a recent paper |Q, Vollmayr-Lee and Luijten (VLL) 
present a Kac-potential approach of nonintegrable inter- 
actions. They consider a d-dimensional classical fluid 
with two-body interactions exhibiting a hard core as 
well as an attractive potential proportional to r~ T with 
< r/d < 1 (logarithmic dependance for r/d = 1) ||. In 
their approach, they also include a Kac-like long-distance 
cutoff R such that no interactions exist for r > R, and 
then discuss the R — > oo limit. They show that the exact 
solution within Boltzmann-Gibbs statistical mechanics is 
possible and that - no surprise (see VLL Ref. [12] and 
references therein) - it exhibits a mean field criticality. 
Moreover, the authors argue that very similar considera- 
tions hold for lattice gases, 0{n) and Potts models. 

Almost as an illustration of the fact that "It is the 
strange privilege of statistical mechanics to stimulate and 
nourish passionate discusions related to its foundations 
[■■■}" B> VLL state "Our findings imply that, contrary 
to some claims, Boltzmann-Gibbs statistics is sufficient 
for a standard description of this class of nonintegrable 
interactions. ", and also that "we show that nonintegrable 
interactions do not require the application of generalized 
q-statistics. ". VLL also inform us that "the main moti- 
vation for this [their] work stems from the considerable 
attention systems with nonintegrable interactions have 
received in the context of "nonextensive thermodynam- 
ics."". 

It is the purpose of the present Comment to argue that 
VLL's (first two) above statements severely misguide the 
reader. Indeed, the interesting VLL discussion, along 
traditional lines, of their specific Kac-like model only 
exhibits that Boltzmann-Gibbs statistical mechanics is 
- as more than one century of brilliant successes guar- 
antees! - necessary for obtaining thermal equilibrium 
properties without needing to do time averages; by no 
means it proves that it is sufficient, as we shall soon clar- 
ify. Neither it proves that wider approaches (such as, for 
instance, nonextensive statistical mechanics, VLL Refs. 



[6,31] and present Ref. or a similar formalism) are not 
required or convenient. The crucial point concerns time, 
a word that nowhere appears in the VLL paper. This is, 
in fact, rather surprising since the key role of t has been 
strongly emphasized in several occasions, for instance in 
VLL Ref. [31] (e.g, Fig. 4 of of VLL Ref. [31] illustrates 
the expectation). For integrable (or "short-range" , as fre- 
quently referred to in the literature focusing the present 
context j|]) two-body interactions in a TV-body classical 
Hamiltonian system, i.e., for r/d > 1, we expect that 
the t — > oo and N — > oo limits are commutable in what 
concerns the equilibrium distribution p{E) , E being the 
total energy level associated with the macroscopic sys- 
tem. More precisely, we expect naturally that 

p(E) = lim lim p(E; N; t) = lim lim p{E; N; t) 
oc exp[-^/fcT] (r/d > 1) (1) 

if the system is in thermal contact with a thermostat at 
temperature T. In contrast, the system is expected to 
behave in a more complex manner for nonintegrable (or 
"long-range") interactions, i.e., for < r/d < 1. In this 
case, no generic reasons seem to exist for the t — > oo 
and N — > oo limits to be commutable, and consistently 
we expect not necessarily equal results. The simplest of 
these results (which is in fact the one to be associated 
with the VLL paper, although therein these two relevant 
limits and their ordering are not mentioned) is, as we 
shall soon further comment, 



lim lim p(E; N; t) oc exp[- 

N—tCXD t—>oo 



(E/N)/(kT/N)} (2) 



where we have introduced N = [N 1 ~ T / d — r/d]/[l — r/d] 
in order to stress the facts that generically 

(i) E is not extensive, i.e., is not proportional to N 
(but is E oc NN instead; more precisely, E is extensive 
if r/d > 1, see || and VLL Refs. [4,5], and it is nonex- 
tensive if < rjd < 1), 

and 

(ii) T needs to be rescaled (a feature which is frequently 
absorbed in the literature by artificially size-rescaling the 
coupling constants of the Hamiltonian) , in order to guar- 
antee nontrivial finite equations of states. Of course, for 
r = 0, we have N = N, which recovers the traditional 
Mean Field scaling. 

But, depending on the initial conditions, which deter- 
mine the time evolution of the system if it is assumed 
isolated, quite different results can be obtained for the 
ordering limt_ >00 limjv-»ooP(-2')-^;t). This fact has been 
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profusely detected and stressed in the related literature 
(see, for instance, VLL Ref. [31], present Refs. j||7]-[ll]] 
and references therein). 

Let us discuss this point further. The impressive suc- 
cess of Boltzmann-Gibbs statistical mechanics relies on 
the fact that, for ubiquitous Hamiltonian systems, time 
averages (whose calculation is typically untractable) can 
be replaced by appropriate ensemble averages (by far less 
harder to calculate) . More specifically, let us consider an 
isolated iV-body system (microcanonical ensemble) for 
which we fix the total energy E (for some systems, other 
constants of motion need to be fixed as well). If we as- 
sume an initial condition for the system and let it evolve 
(following Newtonian mechanics if the system is a clas- 
sical one), after sufficiently long time t (mathematically 
speaking, in the limit t — > oo; practically speaking, for 
times well above the inverse maximum Lyapunov expo- 
nent, whenever it is positive) the distribution of energies 
is given by an expression which is the microcanonical 
equilibrium distribution. From this distribution, we can 
in principle calculate the marginal distribution associ- 
ated to any one among the N particles. This equilibrium 
distribution typically is, as first perceived by Gibbs, a 
power-law (see, for instance, ]^Jl2|). For fixed value of 
E properly scaled with N and N — > oo, this distribution 
approaches the celebrated Boltzmann, exponential factor, 
which corresponds in fact to the canonical ensemble, i.e., 
when the system is in thermal contact with a thermostat 
(which fixes T, of course). A beautiful illustration is ex- 
plicitely worked out in [j/j for the marginal distribution 
corresponding to the momentum of one particle: in the 
N — > oo limit, the celebrated Maxwellian velocity distri- 
bution is neatly recovered. By the way, in 1993 Plastino 
and Plastino remarked |l2|] that the microcanonical equi- 
librium power-law distribution just mentioned precisely 
is the one which emerges from nonextensive statistical 
mechanics where 1 — q plays essentially the role of 1/N. 

What we have just described corresponds clearly to the 
limjv^oo linit^oo p(E; TV; t) ordering. What rele- 
vant modifications are expected to happen in the 
lim t ^oo liniAr^oo p(E; N; t) ordering? Typically none if 
no long-range interactions are involved. But, as men- 
tioned above, the situation is expected to be much more 
subtle in the presence of long-range interactions. To be 
more precise, consistently with the available numerical 
results, we generically expect (as conjectured in Fig 4 
of VLL reference [31], and exhibited in [Q-Q, among 
others) that, after some transient towards equilibration 
starting from certain classes of initial conditions (such 
as waterbag or double waterbag in velocities, for in- 
stance) , the system achieves a quasi- stationary or meta- 
equilibrium state whose duration diverges with N. After 
this state, for finite N, the system typically relaxes to the 
Boltzmann-Gibbs equilibrium (i.e., q — 1). Let us stress 
at this point that, if N ~ 10 23 , during times which could 
be longer than the age of the universe, it is the anoma- 
lous, non Boltzmann-Gibbsian, meta-equilibrium state 
which will be observed and not the standard one (focused 



in the VLL paper). In other words, between the stan- 
dard and such nonstandard macroscopic states, some- 
thing analogous to an activation barrier exists, whose 
height diverges with TV. A typical (very illustrative but 
by no means unique) such system is the classical d- 
dimensional model of localized planar rotators with two- 
body (attractive) interactions which decay as r~ T , and 
whose coupling constant does not (unphysically) depend 
on N . The following anomalies have been (either ana- 
lytically or numerically) observed in the microcanonical 
molecular dynamics approach: 

(i) Above a rescaled critical energy u c = U C /(NN) 
where U denotes the fixed total energy (u c — 0.75 for 
r = 0), the system exhibits a (conveniently scaled) Lya- 
punov spectrum whose largest Lyapunov exponent ap- 
proaches, for increasingly large N, a finite positive value 
for r/d > 1, but vanishes instead for < r/d < 1 (it 
vanishes, in fact, as N~ K , where K,(r/d) seems to be an 
universal function, independent of both d and u, which 
decreases from 1 /3 to zero when r/d increases from zero 
to unity, see VLL reference [35] and present reference 

@); _ 

(ii) Below u c , long standing metaequilibrium states 
are observed everytime the simulation is started from 
single (or even double) waterbag distributions in veloci- 
ties. Specifically, the time evolution of the (properly size- 
scaled) average kinetic energy (also averaged over many 
statistically equivalent sets of initial conditions), exhibits 
]l0| plateaux corresponding to temperatures apprecia- 
bly different from the one associated with the canoni- 
cal Boltzmann-Gibbs distribution. The duration of these 
plateaux seems to diverge with N, while their associated 
height remains finite. This result carries a most impor- 
tant corollary: it strongly suggests that the zeroth prin- 
ciple of thermodynamics can be valid out from the usual 
(q = 1) statistical mechanics Q ; 

(iii) Consistently with point (ii), the one-particle dis- 
tribution of momenta is not the Maxwellian one during 
the entire extent of the plateaux just mentioned. It does 
not change for all times included in the plateaux, it is 
more peaked than the Maxwellian one for low momenta 
(consistently with a nonextensive statistical mechanics 
canonical distribution for q > 1), and exhibits a cut-off 
for high momenta, as generically expected for any mi- 
crocanonical distribution (indeed, since the total energy 
is finite, even if all the energy was momentarily concen- 
trated in the kinetic energy of a single particle, the cor- 
responding momentum could not be infinite); 

(iv) Consistently with points (ii) and (iii), a Levy- 
type anomalous superdiffusion is observed during the 
plateaux, which makes a crossover to normal diffusion 
when the plateau ends, and the usual Boltzmann-Gibbs 
equilibrium is attained. 

It is very clear that no reason at all exists for expecting 
the above list to be exhaustive in describing the anoma- 
lies associated with the meta-stable states. It has been 
given just as an illustration of the important features 
that the VLL paper has overlooked. It is in this sense 
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that we have argued in the beginning of this Comment 
that the VLL paper does not prove, contrary to what 
is therein stated, that Boltzmann-Gibbs is sufficient, it 
only provides one more illustration that it is necessary, as 
probably all mechanical statistical physicists of the world 
are convinced. Do we mean by this that we have herein 
proved that, for classical conservative many-body long- 
range interacting Hamiltonians, nonextensive statistical 
mechanics is also necessary? Certainly not! What we 
believe that we have shown is that another, nonstandard 
formalism is also needed: at the present stage, for the 
particular system we are focusing on, nonextensive sta- 
tistical mechanics is but a candidate (perhaps the correct 
one), on which many scientists are currently working. In 
a kind of desperate last effort, one could easily imagine 
a devil's advocate arguing that Boltzmann-Gibbs statis- 
tical mechanics is the uniquely correct one for describ- 
ing the ultimate equilibrium macroscopic state. On this 
we have no objection at all, it even strongly seems that 
this is a stricly correct statement. However, for long- 
range interacting macroscopic systems initially placed in 
a nonstandard bassin of attraction (in the space of the 
distributions) of the initial conditions, this state might 
be achieved "after the end of the universe" . This feature 
would make its physical utility a very controversial mat- 
ter, of a byzantine style which is not much of our taste. 
Along these lines, the only thing that really matters is 
whether, yes or no, a neat connection can be established 
between the power-law distributions which extremize the 
nonextensive entropic form S q = [1 — X^-P?]/[9 — 1] @] 
and the long standing metaequilibrium states existing be- 
fore the ultimate, Boltzmann-Gibbs equilibrium state is 
attained. To achieve this goal, a tractable way is to per- 
form microcanonical molecular dynamics simulations in 
long-range- interacting systems which contain N » 1 
particles, and then focusing on subsystems of them with 
M >> 1 particles such that N » M » 1. In the 
(M,N/M) — ► (oo, oo) limit the comparison is expected 
to make sense. Work along this possibility is in progress. 

It is certainly relevant to remind at this point that the 
celebrated Boltzmann-Gibbs exponential distribution for 
Hamiltonian systems has been founded in the literature 
in at least four manners, namely the variational entropic 
principle [jl5| , the steepest descent method jl6| , the laws 
of large numbers JT^ ] and the microcanonical counting 
All four have been generalized now (respectively , 
71 and |2l[]), and systematically lead to the power- 
law distributions emerging within nonextensive statisti- 
cal mechanics. 

Although perhaps not totally transparent, the VLL pa- 
per and consistently the present Comment focus essen- 
tially on classical Hamiltonian many-body systems. How- 
ever, the concepts involved in nonextensive statistical 
mechanics have been successfully applied to many other 
systems. These include fully developed turbulence [ p2|j23| 
(Beck has recently developed a quite impressive theory 
with not a single free parameter, which compares remak- 
ably well with the available experimental data), granu- 




lar matter p4j, electron-positron annihilation and other 
high energy systems [|25| , Levy and correlated types of 
anomalous diffusions pq | , low-dimensional nonlinear dy- 
namical systems (^?J , self-organized criticality [^8| , reas- 
sociation in folded proteins |29| , quantum entanglement 
fjof , to quote but a few. For these systems, it is al- 
lowed to think that the question which remains to be 
fully clarified is not whether the formalism works succes- 
fully, but why it does so. Given the wide diversity of the 
systems under focus, it is not yet totally clear what are 
the basic ingredients of the game. It is however already 
acquired that some type of (multi)fractality or hierarchi- 
cal structure is apparently always present (the physical 
mechanisms capable of driving such fractality appear to 
be long-range interactions, strongly nonmarkovian pro- 
cesses, fractal boundary conditions, quantum nonlocality, 
mesoscopic dissipation, and others). This fact tends to 
generate a slow, power-law, mixing in phase and analo- 
gous spaces, as opposed to the usual, exponential mixing. 
The situation might well be, in some cases, more subtle 
than just this (for instance, in the system of long-range 
interacting rotators addressed above, the Lyapunov ex- 
ponents in the N — > oo limit vanish for u > u c but are 
finite for u < u c ; see |3l[]), but it is presently unavoidable 
to think that microscopic mixing issues lay at the heart 
of nonextensive statistical mechanics and thermodynam- 
ics. This possibility goes in fact very well along ideas 
of Krylov ||, Balescu ||, Dorfman @] and others. 
Further studies are certainly needed and welcome. 

It is with pleasure that I acknowledge useful remarks 
by D.H. Zanette, S. Abe and A.K. Rajagopal. 
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To illustrate the kind of subtleties that might be involved, 
let us consider a simple example. The solution of the lin- 
ear ordinary differential equation £ = Ai£ (with £(0) = 1) 
is £ = e Al * , where £ represents the sensitivity to the initial 
conditions and Ai is assumed to be the Lyapunov expo- 
nent of some simple, one-dimensional, nonlinear dynami- 
cal system. If Ai = 0, the appropriate equation might be 
the nonlinear one £ = X q £ q (q £ 7V), whose solution is 
now a power-law, namely £ = [1 + (1 — q)Xqty 1 '^ 1 ~ q ". 
This solution recovers the previous one as the q = 1 
instance. Let us now consider the more general case 
where both linear and nonlinear terms are present. The 
corresponding differential equation can be written as 
£ = Ai£ + (Aq — Ai)£ 9 . The solution is now given by 
a)A 1 t]i/(i-«)_ The pure linear cage 

Ai (Vg) case as well as in 
the q = 1 cases. The pure nonlinear case is recovered 
in the Ai = particular case. An interesting situation 



is recovered in the A, = 
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appears when q < 1 and \ q >> Ai > 0. In this case, 
for t « 1/[(1 — the solution practically coincides 

with the power-law one, and for t >> 1/[(1 — g)Ai the 
solution asymptotically reproduces the exponential be- 
havior. Generally speaking, it is possible to imagine that 
some specific effects could depend basically only on the 
nonlinear term (oc £ 9 ) and practically not on the linear 
term (oc £). If so, the relevant part of would be ba- 
sically related to the power-law behavior, whether Ai is 
zero or not. Considerations of this sort might have rel- 
evance for deeply understanding what happens dynam- 
ically in macroscopic systems like the long-range inter- 
acting rotators one discussed in this Comment, where a 
value u c exists such as at both sides of which unusual 
features are observed: vanishing Lyapunov exponents for 
u > u c , and the existence of long standing metastable 
states coexisting with nonzero Lyapunov exponents for 
u < u c - A situation analogous to this one does occur in 
standard critical phenomena. Indeed, a crucially impor- 
tant quantity (directly related to the correlation length, 
the temperature dependance of which defines in turn the 
critical exponent v) is the two-body space correlation, 
not of the values of the local order parameter at different 
sites, but of the fluctuations, with regard to the average 
order parameter, of the values of the local order parame- 
ter at different sites . The distinction is inexistent for the 
disordered phase (typically at T > T c ), where the aver- 
age order parameter vanishes, but it is most important 
for the ordered phase (typically at T < T c ), where the 
average order parameter is nonzero. 
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